Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN ETeX style file v2.2) 



Radiative pressure feedback by a quasar in a galactic bulge 



A.C. Fabian 1 , A. Celotti 2 and M.C. Erlund 

institute of Astronomy, Madingley Road, Cambridge CB3 0HA 
2 S.I.S.S.A., via Beirut 2-4, 34014, Trieste, Italy 



5 February 2008 



ABSTRACT 

We show that Eddington-limited black hole luminosities can be sufficient to deplete a galaxy 
bulge of gas through radiation pressure, when the ionization state of the gas and the presence 
of dust are properly taken into account. Once feedback starts to be effective it can consistently 
drive all the gas out of the whole galaxy. We estimate the amount by which the effect of 
radiation pressure on dusty gas boosts the mass involved in the Eddington limit and discuss 
the expected column density at which the gas is ejected. An example is shown of the predicted 
observed nuclear spectrum of the system at the end of an early, obscured phase of growth when 
the remaining column density Ah ~ 10 / cm~ 2 , where / is the gas fraction in the bulge. 
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1 INTRODUCTION 

Much observational work over the past decade has shown that the 
mass of the central black hole in a galaxy, Mbh> scales with the 
mass and/or velocity dispersion, a, of the bulge of that galaxy (Ko- 
rmendy & Richstone 1995; Magorrian et al 1998; Gebhardt et al 
2000; Ferrarese et al 2001). A recent correlation (Tremaine et al 
2002) shows that Mbh oc <t 4 holds over at least 3 decades in black 
hole mass. This result suggests that black hole and galaxy growth 
are entwined and introduces the exciting possibility that the growth 
of a central black hole determines the properties of its host galaxy. 

Many models have been produced to explain the Mbh — & 
correlation. The ones relevant here involve a central active galaxy 
influencing the level of gas, and thus star formation, of the host 
galaxy. The black hole may for example grow in mass and power 
until it is capable of ejecting the interstellar medium from the 
galaxy, thus stopping star formaton and determining the total stellar 
mass of the galaxy. Such models either employ an energy argument 
(Silk & Rees 1998; Haehnelt, Natarajan & Rees 1998; Wyithe & 
Loeb 2003) which leads to Mbh oc a or a momentum one using a 
quasar wind (Fabian 1999) or the quasar radiation directly (Fabian 
et al 2002; King 2003; Murray et al 2004), to obtain Mbh oc <t 4 . 
More general heating models have been presented by Granato et 
al (2005) and by Sazanov et al (2005), while Begelman & Nath 
(2005) explored the role of momentum in self-regulating the gas 
density profile. 

The binding energy of the likely interstellar medium of a 
galaxy is less than one per cent of the energy released by the growth 
of its massive central black hole. This energy must of course be 
supplied in order to eject the gas. Sufficient momentum is also es- 
sential, which makes the overall process energetically inefficient. 

Here we revisit the momentum approach using radiation pres- 
sure from a central quasar acting on the dense star-forming inter- 
stellar medium of a young galaxy. The obvious approach is to use 



the Eddington limit, but in its original form it applies to a point 
mass while here we wish to deal with a galaxy which is a distributed 
mass. If the standard Eddington limit cannot be exceeded and is ap- 
plied to both quasar and galaxy, then there is no way that the quasar 
radiation can eject the surrounding galactic medium since the mass 
and thus the Eddington-limiting luminosity required rises with ra- 
dius. One approach (King 2003) is to invoke super-Eddington radi- 
ation levels from the quasar. There is, however, little observational 
evidence for super-Eddington radiation (e.g. Woo & Urry 2002; 
Kollmeier et al 2005) or firm theoretical basis for it (but see Begel- 
man 2002). Here we introduce an effective Eddington limit which 
occurs when radiation acts on lowly ionized and dusty gas. The 
quasar can easily exceed this limit and so drive gas from the galaxy. 
Of particular interest here is the column density of the gas within 
the galactic bulge. 



2 THE EDDINGTON LIMIT 

We now derive an effective Eddington limit for the situation when 
the radiation from a central active galaxy or quasar interacts with 
dusty, partially-ionized gas. The effective interaction cross-section, 
due to photoelectric absorption, dust extinction etc. can then be 
much larger than the Thomson cross-section used to derive the stan- 
dard Eddington limit. A central quasar at the standard Eddington 
limit for its mass, relevant for highly ionized gas in its immediate 
vicinity, can be radiating at, or above, the effective Eddington limit 
for distant matter gravitationally bound by the higher mass of the 
black hole and the host galaxy. 

The Eddington limit arises when the outgoing radiation pres- 
sure, due to electron scattering, from a source of luminosity L bal- 
ances the gravitational attraction due to its mass M: 

i-iTGMm p c 

^Edd = • (1) 
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It is assumed here that radiation pressure is acting on a gas of 
ionized hydrogen around the photon source. G, m p and <jt are the 
gravitational constant, proton mass and Thomson cross-section, re- 
spectively. If the cross-section for interaction between the radiation 
and matter cr\ is larger than for electron scattering, then the rele- 
vant limit which we denote the effective Eddington limit L' Edd is 
proportionally changed. Furthermore, the Eddington limit is often 
derived by considering an isolated electron-proton pair exposed to 
the radiation force. Here we need the effects of radiation on shells 
of matter surrounding the source so we consider the effective limit 
for a column density of gas N. For a gas optically thin to Thomson 
scattering, L' Edd ~ LEddTT/min[Ti, 1], where tt and r, are the 
optical depths for the corresponding cross-sections, r = aN. 

This means that a luminosity which is sub-Eddington for com- 
pletely ionized gas close to a central mass, M, can exceed the mod- 
ified Eddington limit for partially-ionized or neutral gas, for which 
CTi > (Jt, that is denser or further away. Note that a\ is an effective 
cross-section obtained by averaging over the incident spectrum, the 
column density of matter and the state of that matter (ionization 
state, dust content, chemical composition, etc.). It therefore de- 
pends on the spectral shape of the incident radiation. Here we refer 
to the integrated spectrum and thus to an average cross-section. The 
absorbed luminosity correponds to L a ~ Lt\ in the optically thin 
regime and L a ~ L for optically thick gas. We assume i a is radi- 
ated isotropically by the absorber and thus resultant rate of change 
of momentum per unit area, or radiation pressure, is L a /47rr 2 c. 

The amplification factor A owed to the presence of gas not 
fully ionized and dust, can be defined as the ratio of the effective 
radiation pressure Lmin[ri, l]/4ivr 2 c acting outward on a column 
N gas at radius r with respect to that for fully ionized gas. This 
corresponds to 



A = 



(2) 



Alternatively, this can be re-expressed in terms of the mass 
required to gravitationally hold this column density back from ex- 
pulsion, which in terms of A can be written as 



ML 



AM B 



(3) 



The latter expression assumes that the central black hole, of mass 
Mbh, is radiating at its (Thomson scattering) Eddington limit, so 

L = Z/Edd- 



3 THE EDDINGTON BOOST FACTOR 

We have determined the boost factor, A, by finding the luminosity 
absorbed by column density, Nh, after making assumptions about 
the ionization parameter of the gas and incident radiation spectrum. 
For this purpose we use the code CLOUDY 96.01 (Ferland 1993) 
with the AGN spectrum. This essentially assumes a UV blackbody 
plus an X-ray powerlaw spectrum appropriate for a supermassive, 
radiatively-efficent, accreting black hole. 

The boost factor is obtained from the directly absorbed radia- 
tion. We subtract from the input luminosity that which is transmit- 
ted (without any diffuse radiation which is assumed to be isotropic). 
Basically we are using the continuum radiation pressure and assum- 
ing that the absorbed radiation emerges as heat which, on emis- 
sion, interacts either little or with no dynamical consequences on 
the galactic gas. Trapping of radiation is assumed to be negligible. 
The procedure has been stepped over a range of column density Nh 
and ionization parameter £ = L/nr 2 with the results plotted as A 



in Fig. la. Different symbols represent the cases with no dust and 
gas with a Galactic mix of dust. We expect that this covers the con- 
ditions in growing galaxies. If there is vigorous star formation and 
a high metallicity in the gas surrounding growing massive black 
holes then the gas should be very dusty. 

Overall the boosting factor is independent of Nh (see partic- 
ularly the cases without dust), and just given by the ratio of the 
effective optical depth with respect to the Thomson one, until the 
gas becomes optically thick. For larger columns the radiation pres- 
sure does not increase further, while gravity acts on a more massive 
shell, leading to A oc N^ 1 . The normalization of A in the optically 
thin regime and the column corresponding to an effective depth of 
unity clearly depends (inversely and directly, respectively) on the 
ionization state of the gas, parametrized by £. 

As shown below, the effect of dust is to boost the effect of 
radiation pressure by one or more orders of magnitude, which 
- as discussed later - is key to ensure that gas can be depleted 
from a proto-galaxy for Eddington luminosities from the central 
black hole. The transmitted and total spectra for a column density 
Nh = 10 cm~ and ionization parameter £ = 100 is shown 
in Fig. 2 for both the dusty and dust-free cases. It is clear that 
much more radiation is absorbed in the dusty case with the radi- 
ation from the nucleus much reduced at energies above 0.3 eV or 
shorted than 3.6 microns wavelength. This means that the radiation 
pressure from much of the large UV/optical blackbody emission is 
harnessed so causing the threefold or more greater boost factor A, 
when compared with the dust-free case. Since most of the radia- 
tion is absorbed in the dusty case (L a is close to L), Fig. 1 follows 
straightforwardly from equation 2. The precise energy of the re- 
sulting infrared emission bump depends on the dust temperature 
and thus on the radial distribution of the dust. 

The energy of the blackbody emission, assumed to originate 
from an accretion disc around the black hole, depends on the black 
hole mass, shifting to higher energies as the mass reduces. This 
means that the difference between the dusty and dust-free cases 
reduces for lower mass objects (Mbh < 10 s M Q ). The blackbody 
peak is absorbed for all relevant masses for the dusty case, so we 
expect the resulting M — a relation to be robust where the other 
assumptions hold. 



4 EXPULSION OF GAS FROM A GALAXY 

Let us consider now the specific case of a galactic bulge, and in par- 
ticular follow the scenario discussed by Fabian (1999) and Fabian 
et al. (2002). It is assumed that the bulge is isothermal with mass- 
density profile p oc r -2 . The total mass within radius r is given by 



M 



G 



(4) 



of which fraction / is assumed to be in gas, and this has been ac- 
creted onto the black hole, so Mbh = fM. The column density 
exterior to r is then 

/- 2 



N = 



2nGm p r 



(5) 



It should be noted that in Section 2, we treated the gravita- 
tional force as acting on gas concentrated in a thin shell at distance 
r, while the bulge gas is distributed, and this leads to an extra factor 
ln(r max /r) where r max corresponds to the outer boundary of the 
isothermal distribution (Fabian et al 2002). 

As previously mentioned, accretion onto the black hole within 
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Figure 1. a) Top: Boost factor A as a function of column density (expressed 
in terms of neutral hydrogen column -/Vjj) and for varius values of the ion- 
ization parameter log£ from -3 to 3. Circles denote gas with a Galactic 
mix of dust, b) Bottom: Factor l/N^Acr^ = L/L a shown as a function 
of column density. For dusty ionized gas, the factor is close to unity for 
N H > 10 23 cm" 2 and less than 2 for N H > 5 X 10 21 cm -2 . The effect 
of having the dust-to-gas ratio at 0.3 and 0.1 times the Galactic ISM value 
are shown for a column density of 10 23 cm -2 . 



r leaves a column density given by equation (5) beyond. Our treat- 
ment of such a column being concentrated in a shell corresponds to 
the possible scenario where the very same radiation pressure would 
compress the gas into a shell propagating outward. Alternatively, 
one has to consider the extra logarithmic factor: although the value 
of r max is not clearly determined a priori, there has to be a outer 
boundary of the (otherwise diverging) isothermal mass distribution 
oc r~ 2 . Finally, we note that eqn. (5) also arises if the gas within r 



is assumed to be swept up into a shell rather than accreted into the 
black hole. 

The question we want to quantitatively answer here is whether 
radiative feedback from the central accreting black hole is sufficient 
to deplete gas on the galactic scale. 

Gas is expelled from the bulge when the luminosity of the 
accreting black hole L exceeds the modified Eddington luminosity, 
namely A > 1. Or more precisely, gas at radius r is pushed outward 
when M^ dd for the column external to the gas exceeds the mass 
internal to r. Combining eqs. (2), (4), (5) with L = Z/Edd, this 
corresponds to 

4/7TC nG ' 
or from eqs. (3), (4) and (5) 



(6) 



M BH 



/- 4 



TrG 2 m p AN' 



(7) 



The computation presented in Section 3 shows (Fig. lb) that 
ANijt is close to unity over a wide range of column densities from 
10 22 to 10 24 cm -2 . This means that we can replace AN by or -1 
in the above equation so obtaining 



Mb 



/- 4 



nG 2 m ] 



-<xt- 



(8) 



This is the same expression, within a factor of two, as that 
derived on simpler grounds by Fabian (1999), Fabian et al (2002), 
King (2003) and Murray et al (2004). What we have done here is 
to consider a more realistic gas composition and incident spectrum, 
and to quantitatively consider the effect of dust. The result is ob- 
tained principally because most of the relevant incident radiation, 
that in the UV to soft X-ray bands, is absorbed by column densi- 
ties greater than 10 22 cm -2 of dusty gas which is not completely 
ionized. Note that the resulting black hole mass (eqn 8) depends on 
the Eddington fraction of the source as /g dd , so is larger for sources 
operating below the Eddington limit. Kollmeier et al (2005) find in 
their sample of quasars at redshifts z — 0.3 — 2 that most are within 
a factor of ten of the Eddington limit. 

We have assumed that the outer gas is all in a shell at radius r, 
which slightly overestimates the force required. Radiation pressure 
will sweep matter into a dense shell, the column density of which 
will evolve as TV oc r _1 . Provided that it is still within the regime 
where AN is constant, then it will still be driven outward. 

The ionization parameter of matter at radius r when it is about 
to be ejected is 

L2-nGm v 



but at that point L ~ L a 



C /a 2 
, so from (6) 



£ = 8ircm p o~ 



160a 2 2 , 



(9) 



(10) 



where we use a = 100(72 kms" 1 . This means that the line in 
Fig. 1 for £ = 100 is the important one for the general case. We 
envisage that the gas is dusty and so the upper lines in Fig. la are 
most relevant. The input and output spectrum for such a case is 
shown in Fig. 2. 



5 DISCUSSION 

In a very simple model, we envisage that the black hole grows by 
accreting the inner gas. Its mass fM(r) then increases with time, 
its (standard Eddington-limited) luminosity L rises with mass and 
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Figure 2. Input, transmitted and total emergent spectra for a dusty shell with 
column density TVjj = 10 23 cm~ 2 and ionization parameter £ = 100. The 
corresponding spectra for no dust are also shown. The normalization of the 
y-axis is arbitrary. 

the column density N of the gas beyond r decreases until the con- 
dition given by equation (3) is met. The value of A at this point 
is 1// and N fa /JVt. For gas fractions of 10 per cent the final 
column density is then ~ 10 23 cm -2 . The gas is then ejected from 
the host galaxy. 

If the luminosity is sub-Eddington then N ps //Edd-^T- The 
whole mechanism will fail if //Edd is much less than 0.01 at which 
point AN is no longer close to Nt- 

The inner radius of the gas when expulsion occurs is (from eq. 

(5)) 

This is a factor (2/)" 1 times the accretion radius (r a = GM/a 2 ), 
so for a typical value of f — 0.1, r cxp ~ 5r a . The region where 
most of the obscuration occurs is therefore compact. 

Most of the black hole growth will have been by obscured ac- 
cretion, as implied by the observed X-ray background (Fabian & 
Iwasawa 1999; see also Fabian 2004; Brandt & Hasinger 2005; 
Worsley et al 2004; Alexander et al 2005; Civano, Comastro & 
Brusa 2005; Martinez-Sansigre et al 2005) and mid-infrared stud- 
ies with Spitzer (Treister et al 2005). There is already a consid- 
erable population of known luminous AGN with column densities 
A^h ~ 10 23 — 5 x 10 23 cm -2 . In our model the main black hole 
growth phase occurs as the column density reduces to /JVr . 

We do of course require that the black hole continues to be 
fuelled for the gas expulsion time which is several tens of millions 
of years. This could involve a torus and disc around the black hole, 
such as found in many models for AGN (see e.g. Antonucci 1993). 
Anisotropy of the radiation due to such structures will also cause 
the ejection to not be spherically symmetric. We consider that such 
issues are secondary to the basic model outlined above. 

We have assumed that the gas is mostly cold and yet dis- 
tributed throughout the bulge of the host galaxy in a manner similar 



to that of the stars. This requires that the gas is supported in some 
way and we assume that a pervasive hotter phase may be responsi- 
ble. Cold clouds embedded in a hotter medium can drag the hotter 
medium with them when ejected. Why the distribution should be 
p oc r~ 2 is unexplained, although it must roughly occur if the gas 
clouds form into stars which appear to have this distribution. 

The feedback process described here ultimately switches off 
growth of both black hole and host galaxy by radiation from the ac- 
creting black hole interacting with, and expelling, surrounding cold 
gas. Note that for the feedback mechanism to be effective implic- 
itly allows for the cold gas component in the bulge to convert in 
stars. Different feedback processes are expected in massive galax- 
ies, particularly those at the centres of groups and clusters where 
most surrounding gas is hot and ionized and so immune from the 
effects of radiation pressure on neutral and partially ionized effects. 
In such massive objects other processes are required for feedback 
(e.g. Churazov et al 2005) and the slope and normalization of the 
M — a relation may change (Fabian et al 2005). 
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